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Abstract 

Characteristic functions contain complete information about all the moments of a classical distri- 
bution and the same holds for the Fourier transform of the Wigner function: a quantum cfiaracter- 
istic function, or the chord function. However, knowledge of a finite number of moments does not 
allow for accurate determination of the chord function. For pure states this provides the overlap of 
the state with all its possible rigid translations (or displacements). We here present a semiclassical 
approximation of the chord function for large Bohr-quantized states, which is accurate right up 
to a caustic, beyond which the chord function becomes evanescent. It is verified to pick out blind 
spots, which are displacements for zero overlaps. These occur even for translations within a Planck 
area of the origin. We derive a simple approximation for the closest blind spots, depending on the 
Schrodinger covariance matrix, which is verified for Bohr-quantized states. 
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I. INTRODUCTION 

Experiments in quantum optics (see e.g. [ij), atom traps, or other quickly developing 
technologies are rapidly realizing the promiss of a manipulative quantum mechanics. The 
object of developing the field of quantum information has led to a refinement of techniques, 
such that the interference between the wave functions of single atoms, or single modes of 
an optical cavity, can be measured. So far, experiments have mainly been realized for very 
simple states, but theory can anticipate a future where more delicate states, such as the 
superposition of many coherent states or the excited state of an anharmonic oscillator may 
3e made to interfere. Experimental work with Rydberg atoms already point in this direction 

A typical interference experiment superposes two modified copies of the same initial 
state. For instance, in quantum optics, it is easy to achieve the unitary transformation 
that corresponds to a uniform phase space translation (or displacement) of the phase space 
variables x = {p,q). This translated state can then interfere with the original state. In 
general, the unitary translation operator 



exp 



exp 



(1) 



acts on the state to produce the new state = T^\%1)) in strict correspondence to the 
classical translation, x i— t- x + ^, by the chord ^. 22| Thus, given an arbitrary superposition 
of a state and its translation, a\ip) + b\ip^), with |ap + |6p = 1, the probability that this is 
measured to be in the untranslated state is \a + b{ilj\ilj^)\^. 

Evidently, measurements of such probabilities (through repeated preparation) supply 
detailed quantum information concerning these initial states. Better still, the full set of 
possible overlaps defines the complete phase space representation. 
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This is known as the chord function 



quantum optics [ll (or the Weyl function as in 
Wigner function 5|: 



as one of the quantum characteristic functions of 
, which is the Fourier transform of the 
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The latter can be redefined, following Royer [3|, |6|, as 

1 



W{x.) 



\R. 



(4) 



where Ry_, the Fourier transform of the translation operators, corresponds classically to the 
phase space reflection through the reflection centre x, i.e. xq t- 2x — xq. 

Thus, both the Wigner function and the chord function provide complete information 
about states, by telling us how they respond to certain continuous sets of quantum manip- 
ulations. These translation and reflection operators act on Hilbert space in close correspon- 
dence to classical phase space translations and reflections. Hence, their action on excited 
states of anharmonic oscillators should correspond closely to the translations and reflections 
of the Bohr-quantized curves on which they are classically supported. This is well verified 
in the context of the semiclassical theories for the Wigner function [3] and the chord func- 
tion 

a a. 

The Wigner function, W{x.), oscillates with a non-negligible amplitude for all 
reflection centres, x, such that the quantized curve and its reflection intersect transversely, 
as shown in Fig. [TJi. Likewise, all translation chords, ^, for which the translated curve 




FIG. 1: Classical reflection a) and translation b) of a Bohr-quantized curve. Semiclassically, the 
intersection between the original and the transformed curve determines the stationary points and 
the semiclassical phase is given by the area bounded by them (shadow area). 



intersects the quantized curve transversely, are in the region where the chord function, xiO^ 
displays sizeable oscillations. Since the quantized curve is closed, it intersects its reflection 
around x, or its translation by at an even number of points. In both cases, it is the 
area between the pair of curves joining two intersections that determines the phase of the 
oscillations of the Wigner and the chord function respectively, the shadow areas in Fig. [H 
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However, it turns out that in both cases the classical region associated to the quantized 
curve lies on caustics where a simple semiclassical theory breaks down. In the case of the 
Wigner function, this caustic is the locus of reflection centres in the neighbourhood of the 
quantized curve itself, that is, as the centre, x, approaches the curve, pairs of intersections 
with the reflected curve coalesce. This caustic is generic and was already dealt with in 
Berry's original treatment In contrast, the classical region in the space of translation 
chords is the neighbourhood of the origin, whatever the shape of the quantized curve. This 
is highly nongeneric, because all points on the quantized curve intersect the translated curve 
in the limit as ^ — 0, i.e. both curves coincide. Another (nonclassical) caustic, associated 
to the longest chords for which the translated curve still touches the Bohr-quantized curve, 
was already treated semiclassically in jol (along with the corresponding Wigner caustic). 
Thus, the present theory, joining the short chord region to the oscillatory region completes 
the general picture for the chord function of a Bohr-quantized state with a single degree of 
freedom. 

Notwithstanding the difficulty of including the neighbourhood of the origin of chords in a 
semiclassical theory, this region encodes a rich store of information concerning the quantum 
state. On the one hand, the derivatives of the chord function, evaluated at the origin, specify 
all the moments of position and momentum and their products. For a Bohr-quantized state, 
the moments can be identified with classical averages over the corresponding quantized curve, 
which is one justification for considering a Planck area surrounding the chord origin to be a 
classical region. In contrast, one finds points of complete orthogonality between the state and 
its translation within this same neighbourhood, which shows that classical correspondence 
cannot be pushed too far. In the general case where the state has no refiection symmetry, 
orthogonality occurs for isolated points. A theoretical treatment for the pattern of these 
special chords, named blind spots, was presented for arbitrary superpositions of coherent, 
or squeezed states in 10|. One of the objectives here is to extend this analysis to Bohr- 
quantized states. 

Our starting point is a simple integral formula for the chord function that is only valid for 
small chords. This was presented in |8|, but we here provide a fresh rederivation in section 
2 and then go on to show that it leads to classical expressions for the moments. Of course, 
the knowledge of all the moments provides a Taylor series for the chord function, but its 
finite polynomial approximations cannot be joined smoothly to the oscillatory region, which 
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is well described by the standard semiclassical theory in js, 9|. In section 3, we establish the 
presence of blind spots in a neighborhood of the origin for any extended state. We present, 
in section 4, an interpolation that bridges the two regimes of the chord function (small and 
long chords). The full semiclassical theory is compared numerically with the exact result 
in section 5, for an example of a Fock state that is subjected to an unitary transformation 
which breaks its reflection symmetry. Finally, we discuss our results in section 6. 



II. SMALL CHORDS AND MOMENTS OF POSITION AND MOMENTUM 



Consider a general semiclassical WKB state associated to the one- dimensional classical 
manifold defined by / = X) Q, 12], 



d^S,{p,q) 



dqdl 



,jrS,{q,I)+il3, 



(5) 



where S is the generating function of the canonical transformation between the cartesian 
variables and the action-angle variables, namely x = [p, q) (/, ^), the index j enumerates 
the branches of S and /3 is the Maslov correction. The chord function for a WKB state, 
obtained by translating this state and taking its overlap according to ([2]), is given by the 
superposition 

Xu,{i) = Y.^,k{i). (6) 



jk 



where the terms Xjk are given by 



Xjk{$.) = N' dQ 



dqdl 



dqdl 



l-\SjiQ+)-Sk{Q-)-^pQ]+il3j-f^k) 



(7) 



with Q± = Q ± Y- Thus, the oscillatory semiclassical form for each branch of the chord 
function results from the stationary phase evaluation of this integral, to be discussed in 
section 3, but we are here concerned with the limit where the translation is so small that 
the phase between stationary points does not rise above Planck's constant. 

Due to the symplectic invariance of the chord function , we may choose ^ to be parallel 
to the vertical axis, without loss of generality. In this case, = Q- = Q, so that the 
phase difference between the top and bottom branches of the curve is just the curve area as 
a function of Q, which is large and not stationary. Thus, the neglect of these terms leaves 
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only the 'diagonal' terms in (|6]), which are given by 

J2 Xni^) = N'J2j^Q ^ ^^""^ = e-*«''«('')/^ (8) 

Again, making use of symplectic invariance, the right-hand expression can be identified with 
the semiclassical approximation of the chord function for short chords, introduced in Q]: 

X{i) ~ ^gix(e)A^/fi. ^ ^^W)i,-m^p]/f^^ (9) 
Jo 27r Jq 2ti 

This approximation assumes that the classical curve is specified by action-angle variables, 
that is, I{p,q) = X and, conversely, x(6') = {j>{6),q{6)). 

The formula ([9]) holds for any choice of the direction of the small chord ^ Q]. It describes 
the purely classical features of the state, in as much as it is the exact Fourier transform of 
the 'classical approximation' of the Wigner function, Wi^x) = 6[I{p,q) = X]/27r, proposed 
by Berry However, it is more precise to consider this form of the Wigner function as 
a rash extrapolation of the correct form of the small chord approximation to arbitrarily 
large chords. Indeed, it will be here shown that the small chord version encodes quantum 
orthogonalities, as well as classical moments. 

The definition of the chord function ([2]) allows us to calculate the statistical moments of 
p and q in the form of derivatives of the chord function, i.e. explicitly 



(p") = tr p"p = {-ihY ^ 



and (q") = tr q^P = ^ 

^=0 



• (10) 

1=0 



Conversely, if we know all the moments, then we know the chord function, because the 
expansion in a Taylor series of the chord function is 




(11) 



where. 



n + k 

and Pnk are all possible permutations of products of g" and p'^. The equation ( !T2|) corre- 
sponds to the symmetrization of the product q"^p^, being the important feature of the Weyl 
symbols, which guarantees the symplectic invariance of the chord function. 
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According to (fTOj) and (Q, for a given parameterization, the moments are obtained by 
evaluating integrals of the form 

{<f) - / ^ mr and (p") ~ / ^ (13) 

These formulas correspond to the classical expected values of powers of position and mo- 
mentum. They can also be obtained directly from the classical approximation of the Wigner 
function, even though this not so satisfactory in other respects. 



III. NODAL LINES AND BLIND SPOTS 



The chord function ([3]) is the Fourier transform of the real Wigner function and is in 
general complex. Indeed, the fact that it represents a Hermitian operator only implies the 
constraint, x{~^) = x(^)*) where the asterix denotes complex conjugation. On the other 
hand, the cosine and the sine transforms, c(^) and s(^) of the Wigner function are real, i.e. 
the real and imaginary part of x(^)) respectively. In line with the definition ([2]), we may 
construct the hermitian operators 



COS 



^ Ax 

h 



and 



so that c(^) 



and s(^) 



So these are generalizations of the 



potentials cos{kq) and sin(A;g) of cold atoms illuminated by standing waves from lasers 



1^ . Il5| . where k is the wave vector of the laser and q the position coordinate. In the 
case when the state has a centre of symmetry (i.e. there exists a centre, x, such that 
^xlV^) = ±|V^)) then s(^) = 0, so that x(^) = c(^) |8|,|l0|. 



In the general case where there is no reflection symmetry, an intersection of a nodal line 
of c(^) with a nodal line of s(^) defines a blind spot 10|, at which the translated state, 
lip^), becomes orthogonal to \ip) (some examples are illustrated in Figs. [2] and [3]). Because 
x(0) = 1, the origin lies on a nodal line of s(^). In a neighborhood of the origin, the chord 
function may be approximated by 



l-;^^A(x)--l((^Ax)2) + 



(15) 
(16) 



Thus, the nodal line of s(^) crossing the origin is locally parallel to the direction of (x). 



On the other hand, because the origin is a local maximum of the chord function, the 
nodal lines of c(^) for small chords avoid the origin. It follows from f fT6|) that the closest 
nodal line surrounding the origin is given approximately by 

((^ A x)2) = {f)e, + {p')e, - m + m)Uv = = (17) 

if we neglect higher order terms. This positive quadratic form is defined in terms of the 
Schrodinger covariance matrix jl6|, K, which establishes the extent of the state in phase 
space, that is, det K > (27r/i)^ is just the symplectically invariant version of the uncertainty 
principle. The nodal line of c(^) is thus approximated by the ellipse (1171) and the closest 
blind spot lies near the tip of the diameter parallel to (x). It is important to note that the 
present estimate for the pair of closest blind spots depends only on the first and second order 
moments. In the case of a quantized curve treated here, it will verified that the qualitative 
features of the nodal lines are explained by this simple approximation. However, the nodal 
lines of c(^) may show marked influence of the higher order moments, in the case of a 
superposition of coherent states [lO|. 

The highly excited Bohr-quantized states appropriate for semiclassical treatment have a 
covariance matrix that is well described by the classical averages discussed in the previous 
section, such that detK ^ (27r/i)^. Thus, the ellipse f|T7j) lies in the deep interior of a 
neighbourhood of the origin with \/h linear dimensions. This is in line with the discussion 
in 10|: Notwithstanding the delicate quantum nature of blind spots, they can be found in 
the 'classical' neighbourhood of the origin and they are precisely determined by classical 
features. This apparent paradox is resolved by the reciprocal relation between large and 
small scales of pure states in phase space, that follows from the universal invariance of the 
intensity of the chord function for pure states with respect to Fourier transformation {4, 8|: 

lx(0P = ^/rf^lx(^)Pexp|^(£Ar7)| . (18) 

So far, we have only estimated the closest blind spots to the origin. It is hopeless to pursue 
the Taylor expansion f ll6p any further to find further orthogonalities. On the other hand, the 
real and imaginary parts of the short chord approximation 1^ have many nodal lines in the 
region where it holds, in the case of a big Bohr-quantized state. In the following section, an 
interpolation formula is presented that allows for a uniform description, which is also valid 
in the outer oscillatory region, where the chord function is evaluated by stationary phase. 
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IV. JOINING THE LONG AND SHORT CHORD REGIMES 



In section 2 we rederived the 'classical' approximation ([9]), which holds in the neighbor- 
hood of the origin, but not in the region beyond. Otherwise, the long- chord regime is well 
described in terms of Airy functions 9] , describing an oscillatory region in a ring surrounding 
the origin, through an outer caustic and on to an asymptotic evanescent regime as |^| — t- cx3 
. This Airy function results from the uniform approximation based on the stationary points 
of the exponent in ([7]). As mentioned in the Introduction, the stationary points are geomet- 
rically identified by the intersections between the supporting manifold and its translation 
{see fig. H]). 

There are two basic reasons why the standard uniform approximation technique, involving 
a transformation of the integral (|9]) into a simpler one (cf. 17|, |l8|), cannot be employed near 
the origin. One is that there are an infinite number of stationary points of the exponent ([7]) 
at ^ = 0, whence the origin is a non-generic caustic (in the sense of the Thom's classification 
theorem jl^ Iisl). On top of this, the large parameter condition in the exponent, essential 



for asymptotic expansions of integral [18[, is not fulfilled for the case of small values of |^|, 
where the behavior of the integral is given by dH]), because \$\/h ~ 1. 

In order to describe the transition regime (between short and long chords), we propose 
the following (semiclassical) expression for the chord function: 

xscio = xsio - spixsm + sp[x^m. (19) 

Here Xs is the integral (|9]), Xw is the full semiclassical integral for the chord function, ([6]) 
and ([7]), and SP[-] denotes the approximation of [■] by stationary phase. Notice that Xs{$) 
was derived as a short chord approximation to Xw{$,), but, since $, A x(6') in the exponent of 
([9]) is a nonlinear function of the integration variable, 6, this integral can also be evaluated 
by stationary phase. Indeed, the pair of integrals in ( |T9l) that are evaluated by stationary 
phase must cancel in the neighborhood of the origin. On the other hand, the middle term, 
SP[xs], cancels Xs for large values of where stationary phase evaluation is valid. 

The stationary phase method assumes that the integral is dominated by points where 
the phase is stationary. In the case of the chord function, such stationary points have a 
geometrical interpretation. Namely, each stationary point defines pairs of points x± on the 
curve. The point x_ is the intersection of the classical curve with its translation by the 
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vector — ^, whereas x_|_ = x_ + ^. This fitting of the chord into the curve is called a chord 
realization. The stationary points of the chord function are the g-coordinates of the centers 
Xj of such realizations. In fig. ^ we show that each chord has two realizations in a convex 
closed curve. 

The amplitude in the above semiclassical approximation may be expressed in terms of the 
canonical action variable, /(x). Specifically, the amplitude is given by |{J+, J~}|~5j where 
is the action variable at the tips of the realization of ^. Finally, the phase in this 
approximation is determined by the area Ay^-, between the realization of ^ and the curve 
plus the product Xj A An additional phase, Uj, is given by the sign of the expression 

^[5(x,+|/2)-5(x,-|/2)]. (20) 

Therefore the stationary phase evaluation for Xw js], is 

SP\X.U) = - E ^^^^■''^ ^ ^ ^ ^ . (21) 

' l{/;.-r7}l' 

where j enumerates the stationary points Xj and 7^ is the difference between the Maslov 
corrections in ([7]). 

On the other hand, for Xs we have 



27r ^ V27rh 



^x(^,(0)A^ + ^sign (x"(^,(0)AO 



(22) 



where i enumerates the stationary points for Xs, that is, the points where the vector ^ is 
tangent to the curve. Thus, inserting the explicit stationary phase evaluations (pT]) and (!22|l 
into f|T9l) leads to a general approximation, in which the only integral left to be evaluated 
numerically was already present in the short chord approximation (Q. 



V. NONLINEAR EVOLUTION OF A FOCK STATE 

In order to test the general approximation f|T9l) . we analize the chord function of a one- 
parameter set of states, evolving under the action of a simple cubic Hamiltonian, depending 
only on momenta: 

H{p) = asp^ + a2p'^ + aip + ao. (23) 
The classical evolution is then determined by 

p{t) = p{0) and q{t) = g(0) + (Sagp^ + aaP + ai)t. (24) 
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For a fixed time parameter t > 0, the action function is I{p,q,t) = {q — [3a3p^ + 2a2P + 
p"^, so tliat tlie classical curve supporting tlie evolved state corresponds to tlie level 
curve I{p, g, t) = X. 

For t = 0, we choose a Fock state, \n) ^ which is an excited state of the harmonic oscillator 
of frequency w. Choosing = 1, Fock states are supported by a circular manifold, so they 
are symmetric with respect to reflections. By choosing__the center of the supporting circle as 
the origin, we obtain the (exact) real chord function 



where i^n(') is the rath Laguerre polynomial 20|. Since Xn is real, it has full nodal lines 



h 



(25) 



corresponding to circles, and its radii are h times to the roots of the nth Laguerre polynomial, 
as shown in fig. |2l The above classical evolution breaks the original central symmetry. 




FIG. 2: (Left) Chord function for the 5th Fock state. It has reflection symmetry and its five nodal 
lines correspond to circles. (Right) Radial cut of X5{^)- Here h = 0.1 



because of the cubic term in the Hamiltonian. For these states, the approximation for the 
first nodal line, eq. (fT7|) . provides a circle of radius ^h/{n+ |). This result differs to the 
exact with ~ 18% of accuracy; so as previously mentioned, the approximation ( IT71) just gives 
a qualitative estimation for the first nodal line. 

On the other hand, the exact evolving chord function, obtained from ([2]), is 



x{i.t) = m)\T^^m)) 



dpC(p+)e-^(*[^(^'+)-^(^-)l-^«'Vn(p-), 



(26) 
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FIG. 3: The chord function for the evolved Fock state e~'^^^^^\n = 5). At this resolution, the 
difference between the exact and the semiclassical approximation is indiscernible. The black lines 
in the real a) and imaginary h) parts are the levels curves for x(^) = 0- c) is the logarithmic 
intensity and d) is the phase. The points in a) show the location of the 'blind spots'. Here /i = 0.1, 
t = 0.1, a2 = Oil = 1. 



where p± = p ± and 



Mp) = ^7= hJ^] e-^ (27) 



is the Fock state n in the p-representation and -&„(■) is the nth Hermite polynomial [20|. 

For any t > 0, the reflection symmetry is broken, thus the chord function has no more 
a global parity, which implies that the imaginary part is not null. As shown in the figures 
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[31 the expression (fT9i) describes appropriately the behavior of the chord function for short 
and long chords. Cleary SP[xw] is not valid when ^ tends to a diameter, but the resulting 
singularity is avoided by using the uniform approximation We can recognize the real (fig 
[3^) and imaginary (fig [2b) parts by their even and odd parity, respectively. The intersections 
of nodal lines for both the real and the imaginary parts correspond to the blind spots, i.e. 
the zeroes of the intensity (local maxima in fig [St) and singuralities in the phase (fig[3li). 
The pair of closest blind spots, which were approximately specified by the first and second 
order moments in section 2, result from the intersection of the smallest closed curve in fig [3^ 
with the straight line through the origin in fig [3)3. It is curious that the other intersections 
also seems to imply radial straight lines in this example. 

The detailed comparison of the intensities for our semiclassical approximation of the chord 
function f[T9|) with the exact result is shown in fig[ll The particular radial straight line chosen 
exhibits a sequence of blind spots. We observe a caustic (the diameter singularity) at the 
left edge of the figure. It arises in the stationary phase approximation, but this spurious 
divergence is corrected by the uniform approximation in j^. 



1\ 

t 




FIG. 4: Comparision between the exact and the semiclassical intensities of the chord function for 
a evolved Fock state e~^^^^^^'\5), along the line = 0.8172^^. We can observe a sequence of blind 
spots. Here h = 0.1, t = 0.1. 
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VI. DISCUSSION 



The chord function portrays a pure state by exhibiting its overlap with all its possible 
translations. There are two cases where such a state has a clear classical correspondence, 
the superposition of well separated coherent or squeezed states treated in {q] and the Bohr- 
quantized states analyzed here. Contrary to naive considerations, in neither case is the decay 
in the square modulus of the overlap smooth and classical like, within a Planck area of the 
chord origin. We have here shown that blind spots, denoting zero overlap, arise deep within 
this classically small neighbourhood. This feature depends basically on the Schrodinger 
covariance matrix: The greater its determinant, the closer to the origin will the blind spots 
lie. 

Knowledge of a finite number of moments determines the chord function of a Bohr- 
quantized state near the origin, but it is insufficient to follow through the complex oscilla- 
tions, punctuated by a complex pattern of blind spots, up to the outer limit of an evanescent 
region. We have presented a new semiclassical approximation for the chord function of Bohr- 
quantized states and verified that it is accurate, right up to the outer caustic, which was 
previously treated in {q]. In the absence of refiection symmetry, such as present in a Fock 
state, the chord function is fully complex, which leads to richer structure than that displayed 
by the corresponding Wigner function. 

If the system is in contact with an uncontrolled environment, the state will not remain 
pure, so it must be described by the density operator p and the definition ([2]) becomes 
x(0 = tr pT-^. Though this chord function still supplies a complete description of the 
state, the overlap with its translation is now given by 

m ^ tr [p{%pfl)] = ^J Ixm^e-'^'^'^'drj. (28) 

For small displacements, ^ — >■ 0, the correlation has a maximal value, since C(0) = tr = 1. 
As the displacement increases, we attain an oscillatory regime, where the stationary phase 
approximation takes account. This behavior changes after the diameter caustic is reached, 
followed by an evanescent region. For the pure states, this correlation is invariant under 



Fourier transform, since it reduces to C(^) ~ | ('?/;|T^|'?/') p = 
As in the case of superpositions of Gaussians wavepackets 



according to (ITS]) . 



10|], we expect that blind 



spots do survive in the chord function for a markovian quantum evolution appropriate to 
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an open system. However they should disappear from C(^) much more quickly than the 
negative regions of the Wigner function. Blind spots correspond to sharp indentations on 
a background of maximal correlations which makes them measurable. As shown in j^, 
they take their place as very sensitive indicators of the full quantum coherence alongside 



the zeroes of the Husimi function 
regions 



2l[ |. The latter often has its zeros in shallow evanescent 



21| . where they are tricky to distinguish. 



Acknowledgments 

We thank partial financial support by CNPq, FAPERJ, INCT-IQ Informagao Quantica 
(brazihan agencies) and CAPES/COFECUB. 

References 



[1] Leonhardt U 1997 Measuring the quantum state of light (Cambridge: Cambridge Univ. Press) 
[2] Haroche S and Raimond J M 2006 Explore the quantum: atoms, cavities and photons (Oxford 
Univ. Press) 

[3] Ozorio de Almeida A M 1998 Phys. Rep. 295 265 

[4] Chomitasis S and Vourdas A 1998 Phys. Rev. A 58 848 - 855 

[5] Wigner E P 1932 Phys Rev 40 749-759 

[6] Royer A 1977 Phys Rev A 15 449 

[7] Berry M V 1997 Phil. Trans. R. Soc. 287, 237 

[8] Ozorio de Almeida A M, Vallejos R and Saraceno M 2005 J. Phys. A: Math. Gen 38 1473-1490 
[9] Zambrano E and Ozorio de Almeida A M 2008 Nonlinearity 21 783-802 

[10] Zambrano E and Ozorio de Almeida A M 2009 New J Phys 11 113044 (14pp) 

[11] Van Vleck J H 1928 Proc. Math. Acad. Sci. U.S.A 14, 178 - 188 

[12] Maslov V P and Fedoriuk M V 1981 Semiclassical Approximation in Quantum Mechan- 
ics (Reidel, Dordrecht) (translated from original russian edition, 1965). 

[13] Ozorio de Almeida 1988 Hamiltonian systems: Chaos and quantization (Cambridge: Cam- 
bridge Univ. Press) 

15 



[14] Harper P G 1955 Proc Phys Soc A 68 874878 
[15] Aubry S and Andre G 1980 Ann Israel Phys Soc 3 133 
[16] Schrodinger E 1930 Proc Pruss Soc Acad Sci 19 296 
[17] Berry M V 1976 Adv. in Phys. 25 1, 1-26 

[18] Dingle R B 1973 Asymptotic expansions: Their derivation and interpretation (Academic Press: 
London) 

[19] Thom R 1975 Strcutural stability and morphogenesis (Reading, MA: Benjamin) 

[20] Abramowitz M e Stegun I 1964 Handbook of Mathematical Functions (New York: Dover) 

[21] Toscano F and Ozorio de Almeida A M 1999 J. Phys. A 32, 6321-6346. 

[22] In the optical context is usually referred to as the displacement operator and is expressed in 
terms of creation and annihilation operators for the harmonic oscillator. This is inconvenient 
for semiclassical analysis. 

[23] An alternative interpretation is to consider c(^) and s(^) to be the chord functions for appro- 
priately symmetrized states [lo| . 



16 



